Suppression of Ecological Competition by Apex Predator 
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In the framework of Lotka-Volterra dynamics with evolutionary parameter variation, it is shown 
that a system of two competing species which is evolutionarily unstable, if left to themselves, is 
stabilized by a commmon predator preying on both of them. Game-theoretic implications of the 
results are also discussed. 
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From the dominant plants in forest vegetation to the 
wild beasts in savanna, an often encountered ecological 
paradox exists in the form of stabilizing influence of the 
top predator. While two species in direct competition 
rarely form a stable ecosystem, they often coexist under 
the dominance of a common predator .1] . The most illus- 
trative example is found in the trophic pyramids, where 
the apex predator, the most savage aggressor of all, seem 
to act as the key guarantor of the stability of the whole 
system Q. For species with intellectual capacity, this 
fact might be explained as a result of conscious effort of 
enlightened self-interest. But the peacekeeping function 
of the apex predator is so prevalent throughout ecosys- 
tems, that the existence of a simple and universal dy- 
namics should be suspected. 
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The purpose of this paper is to understand the struc- 
ture and stability of ecosystems composed of competing 
species in the framework of evolutionary population dy- 
namics 0, ■ The tool we employ is the Lotka-Volterra 
equation with adiabatic parameter variation In this 
approach, the ecological dynamics is determined by the 
time variation of the variables representing the popula- 
tion of the species, while the adiabatic parameter vari- 
ation represents the behavioral evolution of the species. 
The viability of a species in this framework is judged 
both by the short-time ecological stability of the orbit 
and also by the long-term evolutionary stability of the 
shifting parameters. 

We focus specifically on a system that consists of two 
self-sustaining but competing species and an apex preda- 
tor who preys over both competitors. We show that the 
system evolves towards an evolutionarily stable config- 
uration in which the warring preys are tamed into the 
peaceful coexistence. This is in contrast to the case of 
two competitors left to themselves, in which there are no 
evolutionarily stable solution for coexistence, and "arms 
race" drives one of the competitors into eventual extinc- 
tion. We also show that our results can be interpreted in 
a game-theoretic language as the apex predator turning 



FIG. 1: A symbolic diagram showing the pecking order of the 
three species system described by (1). The arrows represent 
the aggression and predation with the specified intensities. 



the prisoner's dilemma between two competitors into a 
collaborative game. 

Suppose there are two populations of competing 
species x\ and x 2 which are preyed upon by a common 
predator population y (Fig. 0). We describe the eco- 
logical dynamics of this system by the standard Lotka- 
Volterra equation 



X\ = b\X\ — a\x\ — a 2 x\x 2 — pixiy, 
x 2 = b 2 x 2 - a 2 x 2 - o\X\Xi - p 2 x 2 y, 
y = -dy + fpixiy + fp 2 x 2 y. 



(1) 
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Here, 61, b 2 are the reproduction rates for species x\, x 2 , 
and a±, a 2l the environmental limitation factor to their 
growth. The coefficient d is the decay rate for the preda- 
tor y, and /, the efficiency of its predation. In the last 
equation, a term proportional to y 2 could be added for 
consistency with other equations, but this can be shown 
to introduce simply a technical complication without af- 
fecting the main line of our arguments. Also, specifying 
separate predation efficiencies for X\ and x 2 makes no es- 
sential difference, because the result can be turned into 
the original form (1) with rescaling of variables. The pa- 
rameters pi and p 2 are the aggression intensities of the 
apex predator y towards X\ and x 2l respectively. Sim- 
ilarly, <ji and 02 are the aggression intensities of x\ to 
x 2 and of x 2 to x\ , respectively. We assume all parame- 
ters to be positive real numbers. At this stage, we treat 
all of them as fixed numbers, making no distinction be- 
tween the Roman denominated "environmental" parame- 
ters and Greek denominated "behavioral" parameters. A 
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nontrivial fixed point Xi (t) 
with i = 1, 2 is given by 
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y (t) = y {x t = y 

d(a 2 pi - cr 2 p 2 ) ~ f{b 2 Px - Wpi)P2 

V d\p\ - (<7l + <T 2 )piP2 

C1P1) + fihpi - b 1 p 2 )pi 



a 2 p\ 



0) 



(2) 



Xi,X 2 



d(a 1 p 2 



a.2Pl 



Y = — — 



-a x p\ 
aia 2 - 



- [<7i + a 2 )PiP2 
a\a 2 



a 2 p\ + a x p\ - {a x + a 2 )pxp 2 



| (q 2 bi - b 2 cT2)pi + (ai&2 - bxax)p2 
a.2p\ + a x p\ - (cri + a 2 )pxP2 

The stability of the fixed point is determined by the be- 
havior of the linearized map 



-axXx -er 2 Xi 
M = I -axX 2 -a 2 X 2 -p 2 X 2 
fpiY f P2 Y 



(3) 



Namely, the fixed point is stable when real part of all the 
eigenvalues A of M determined by 



\XI-M\ =0 



(4) 



is negative. 

When the fixed point is of stable, attracting sort, 
neighboring orbits form an absorbing spiral in phase 
space. We now assume that evolutionary pressure of se- 
lection and adaptation are at work. We can then regard 
aggression intensities px , P2 , o\ and a 2 as evolutionarily 
adjustable parameters which evolve along the path that 
simultaneously increase the functions X 2 \a 2 \ and 

y[pi,/92] until they reach the optimal values. There are 
several indirect pieces of evidence supporting the exis- 
tence of this type of adiabatic evolution among real-life 
ecosystems 0, Q ■ It is convenient to start with the max- 
imization condition for the apex predator dY/ dp\\ p * = 



and dY /dp 



. We then have the relations 
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d 2aia 2 bx 
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bia+ax + aib 2 a- 



P2 = 7 



a 2 b\ - 
2a\a 2 b 2 



a\b\ - &1&20 + 
- b 2 a + a 2 - a 2 bx<T- 



(5) 



a 2 bf + axb 2 - bxb 2 a+ 
These conditions give the expressions 
2a 2 b\ - b 2 a + 



XI 



x 2 



Y* 



4axa 2 — o\ 
2axb 2 — bxaj 



(6) 



4axa 2 - ai 



d 



+ a\b\ — b 1 b 2 a + 
Aaia 2 - a\ 



The quantities — X\ and — X 2 as functions of ax and a 2 
act as the "potential surface" for the variation of a\ and 
a 2 . In (5) and (6), the notation a± = ax ± 172 is used. 
With the definitions a = yjaxa 2 and /3 = \J a 2 /ax • bx/b 2 , 



X<0 



X>0 



2 Oi+ai 



ap/2 a/p ap 



FIG. 2: Fixed point coordinates X* and X% as functions 
of <ti + (72. The parameters are ai = 1, 61 = 1, Gt2 = 1.2, 
62 = 0.9 and d = 2. Solid line represents X* and the dahed 
Xf. The fixed point is stable in the region below a//3 but 
unstable above a/3. The region in between is unphysical. 




FIG. 3: Fixed point coordinates X* and XJ as functions of 
o\ and <72. The parameters are a\ = 1, 61 = 1, = 1.2, 
62 = 0.9 and d — 2. In the region ai + (72 < a C r ~ a/3/2 
= (22&1/62, both cri and (72 have to be decreased to make X\ 
and X2 larger. 



valid parameter range for Xi and Y being positive and 
stable (3?A < 0) is given by 



ax+ a 2 < min {a/3, a/ /3}. 



(7) 



That stability requirement is satisfied can be checked by 
the fact that all the coefficients of the third order polyno- 
mial equation (4) are of same sign within this parameter 
range. 

The evolution of ax and 02 depends on their starting 
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FIG. 4: Aggression intensities pi and p\ as functions of cri 
and <J2- The parameters are ai = 1, b\ = 1, 0,2 = 1.2, &2 = 0.9, 
d = 2 and / = 0.7. 



values. With a straightforward calculation, we obtain 



dXt dXi, 
< < 



do\ da 2 



: (Tl+<T2<0' C , 
< : 0\ + (72 > (T c 



(8) 



within the range of (7). The critical aggression intensity 
a cr is given by 



: {a/3/2, a/2/3}. 



(9) 



If the sum cti +02 is bellow cr cr , both ci and 02 will move 
toward cri = (T2 = 0. Namely, two competing species shall 
settle for a peaceful coexistence as the common preys 
of a predator y. On the other hand, if the sum starts 
above critical value, o\ and a 2 will increase until one of 
the competing species is extinct at that critical value. 
The situation becomes immediately clear with a glance 
at numerical example depicted in Fig. [21 and Fig. [3] 

A crucial point is that the master y acts as a punisher, 
according to (5), that inhibits the increase of a\ and a 2 . 
Fig. 0] serves as a graphical illustration of this effect; In- 
creasing <ii will induce an increase of pi that incur the 
damage upon x\. We stress that no special mechanism is 
assumed for y to police the system in the outset, yet the 
dynamics seems to explain our common sense observation 
about apex predators. 

An intriguing fact is that the critical value a cr is in- 
versely proportional to the natural population of one of 
the prey species, b\j a\ or 62/02- This means that the co- 
existence of competing species under common predator 
becomes a less likelier outcome for a system with richer 
resources. This seems to give a partial explanation to the 
experimentally observed decrease of species at the base 
levels of trophic pyramids [||. 

We next consider the case where the predator leaves 
the scene, namely y = (Fig. [SJ. By setting pi = p 2 
= 0, we obtain, in place of (6), 



XI = 



X* = 



a 2 bi — a 2 b 2 



a\a 2 
aib 2 



ai<j 2 



(10) 



aia 2 — <Ji<t 2 

The linearized map now takes a two-by-two matrix form 



M = 



-a 1 X 1 
-<JiX 2 



-a 2 X x 

-a 2 X 2 



(11) 



<J 2 



FIG. 5: A symbolic diagram showing the two competing 
species described by (1) with pi — p2 = 0. 




FIG. 6: Fixed point coordinates XI and X| for the case 
without the apex predator as functions of <ti and 02- X\ will 
increase o\ and will increase 02 to achieve local advantage, 
which will eventually doom X| to extinction. The parameters 
are a\ = 1, bi = 1, 02 = 1.2, and &2 = 0.9. 



in place of (3). The straightforward calculation gives 
the condition for X± and X 2 to be a viable fixed point, 
namely, X±, X 2 > 0, 5RA < 0, in terms of the allowed 
region for the aggression intensity as 



<7i < min{a,a/3}, a 2 < min{a, a//3}. 



(12) 



However, within this region, we can easily check the re- 
lation 



>0, 



da 2 



> 0. 



(13) 



Therefore, in this case, both o\ or <j 2 shall eventually 
be increased beyond the range (12), and there is no evo- 
lutionarily stable coexisting solutions for two competing 
species. Namely, in the absence of the common master, 
one of the competing species is always driven to extinc- 
tion by arms race of increasing o\ and <j 2 . An example 
of this case is illustrated in Fig. 

In order to fully understand the generic shape of the 
ecosystem, we would have to generalize our arguments 
to more trophic levels than two, and also more species 
than two within a single trophic level. To this end, this 
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TABLE I: The game table X* [<ri, 02] for x\ discretized at 
<7i, o 2 = a/3 (low aggressiveness, "dove") and 2a/3 (high ag- 
gressiveness, "hawk"). Left hand side is the table for the case 
with apex predator (6), and the right hand side, without apex 
predator (10). The game table X|[oi,02] for X2 is obtained 
by transposition of raw and column. The Nash equilibrium is 
indicated with boldface. 

With Apex Predator No Apex Predator 



CT1 \ CT2 


a/3 
dove. 


2a/3 
hawk 


\ "2 


a/3 
dove 


a/3 


3b 


36 


a/3 


66 


dove 


8a 


9a 


dove 


8a 


2a/3 


3b 


36 


2a/3 


66 


hawk 


9a 




hawk 


7a 



work, in combination to our previous work, is but a mod- 
est start. More involved and sophisticated approaches of 
both numerical and analytical nature may have to be 
called for 0, 0, 0| • In the current work, no precise spec- 
ification for the evolutionary dynamics of behavioral pa- 
rameters has been required. While we stress that this is 
an advantage, construction and analysis of more detailed 
models with such specifications are attractive possibili- 
ties. 

Final ly, s ome remarks in the broader context of game 
theory [lfj are in order. There is an obvious game- 
theoretic interpretation of the results (6) and (10). For 



the sake of simplicity, let us set a% — a 2 = a and 
bi = b 2 = b. We now regard (i — 1, 2) as the con- 
trol parameter of the strategy of population a;, for the 
game played between x\ and xi whose payoff tables are 
given by X*[a\, 02] and X£ [<7i, a 2 \. To facilitate the un- 
derstanding, the game tables for two discretized points 
for o"i and 172 are tabulated in Table I. For the case with- 
out the common predator, (10), the game is a continuous 
strategy version of prisoner's dilemma 0, whose 
outcome is o\ — > a, 02 ~~ * a which leads to the extinction 
of either x\ or x 2 - With the introduction of the apex 
predator, the game table is turned into one of collabo- 
rative game, whose outcome is the coexistence a\ — > 
and o"2 — > 0. Note that the game table is symmetrized 
under the presence of the apex predator; X 2 [a\,a 2 ] = 
X*[a 2 , o"i]. This could be interpreted as the sign of altru- 
istic behavior . Advantage of having the apex preda- 
tor as a "law enforcer" is evident, and the loss of half of 
the populations to the predation would be an acceptable 
tradeoff. 

In summary, we have established, for Lotka-Volterra 
systems with evolutionary parameter variation, that two 
competing species are evolutionarily unstable, but can be 
stabilized by the introduction of an apex predator. We 
hope this to be a start for systematic understanding of 
stable ecosystems. 

The authors wish to express gratitudes to Professors 
Kazuo Takayanagi, Toshiya Kawai and David Greene for 
helpful discussions and useful comments. 
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